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Abstract 

The channeling of the ion recoiling after a collision with a WIMP changes the ionization signal in direct 
detection experiments, producing a larger scintillation or ionization signal than otherwise expected. We give 
estimates of the fraction of channeled recoiling ions in solid Xe, Ar and Ne crystals using analytic models 
produced since the 1960's and 70's to describe channeling and blocking effects. 
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1. Introduction 

Channeling and blocking effects in crystals refer to the orientation dependence of charged ion penetration 
in crystals. In the "channeling effect" ions moving along symmetry axes and planes in a crystal, suffer a 
series of small-angle scatterings that maintain them in the open "channels" in between the rows or planes 
of lattice atoms and thus penetrate much further into the crystal than in other directions and give all their 
energy into electrons. Channeled incident ions do not get close to lattice sites, where they would be deflected 
at large angles. The "blocking effect" consists in a reduction of the flux of ions originating in lattice sites 
along symmetry axes and planes, due to large-angle scattering with the atoms immediately in front of the 
originating lattice site. 

In the context of direct dark matter detection channeling occurs when the nucleus that recoils after being 
hit by a WIMP (Weakly Interacting Massive Particle) moves off in a direction close to a symmetry axis 
or symmetry plane of the crystal. Non-channeled ions loose most of their energy into lattice atoms with 
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which they collide. Thus channeled ions produce more scintillation and ionization than they would produce 
otherwise, since they give all their energy to electrons. 

Ion channeling in crystals has recently received a large amount of attention in the interpretation of direct 
dark matter detection data. The potential importance of this effect was first pointed out for Nal (Tl) by 
Drobyshevski 1] and then by the DAMA [2j collaboration, which estimated the fraction of channeled recoils 
to be very large, close to 1, for low recoiling energies in the keV range. This estimate of the channeling effect 
leads to considerable shifts in the cross section versus mass regions of acceptable WIMP models in agreement 
with the DAMA/LIBRA [3J annual modulation data towards lower WIMP masses. As a consequence, the 
comparison between the DAMA results and the null results of other experiments was affected too (see e.g. 
Ref. ^| and references therein). 

In Ref. [5] we showed that the blocking effect, not taken into account in the DAMA evaluation, is impor- 
tant and thus the fraction of recoiling ions that are channeled is smaller than in the DAMA collaboration 
estimate. The nuclei ejected from their lattice sites by WIMP collisions are initially part of a crystal row 
or plane, thus, as argued originally by Lindhard 6] (through what he called the "Rule of Reversibility"), in 
a perfect lattice and in the absence of energy-loss processes, the probability that a particle starting from a 
lattice site is channeled would be zero. However, any departure of the actual lattice from a perfect lattice, 
e.g. due to vibrations of the lattice atoms which are always present, violates the conditions of this argument 
and allows for some of the recoiling lattice nuclei to be channeled. In Ref. Q we showed that the channeling 
fractions for Nal (Tl) crystals would never be larger than a few percent. This implies that channeling does 
not change the regions in cross section versus mass of acceptable WIMP models in agreement with the 
DAMA/LIBRA annual modulation data at less than 7a jj. 

Channeling in dark matter direct detection might in principle also lead to a new background free dark 
matter signature. This was pointed out initially in 2008 and later in 2010, by Avignone, Creswick, and 
Nussinov |8| who suggested that a new type of daily modulation due to channeling could occur in Nal (and 
other) crystals measuring ionization or scintillation. Such a modulation of the rate due to channeling is 
expected to occur at some level because the "WIMP wind" arrives to Earth on average from a particular 
direction fixed to the galaxy. Thus, Earth's daily rotation changes the direction of the WIMP wind with 
respect to the crystal axes and planes, changing the amount of recoiling ions that are channeled vs non- 
channeled, and therefore the amount of energy visible via scintillation or ionization. If this daily modulation 
could be measured, it would be a signature of dark matter without background since no other effect could 
depend on the orientation of the detector with respect to the WIMP wind. Following this idea, in Ref. Q we 
provided a full computation of the effect and using our estimates of the channeling fractions in Nal obtained 
upper bounds to the expected amplitude of daily modulation due to channeling in Nal crystals. We found 
large daily modulation amplitudes of the signal rate, even of the order of 10% for some WIMP candidates. 
However, even the largest amplitudes would not be observable in the 13 years of combined DAMA/Nal 
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and DAM A/LIBRA 3] data because of their large background. In the hope that the daily modulation due 
to channeling could be observable in future dark matter experiments, we have provided estimates of the 

. n n 

channeling fractions of recoiling ions also in Ge, Si [1 01 ] . and Csl [ll| crystals, and in the present paper we 
study noble gas crystals. 

The level of background in future experiments is a crucial element to determine if the daily modulation 
is observable and the advantage of u sing solidified noble gas detectors is the possibility of achieving a low 



background level for WIMP searches |12| . Xenon and Neon have no long life radioisotopes and thus contain 
no intrinsic background source of radiation. Argon has the drawback that it contains 39 Ar beta source which 
induces electron recoil signature in the detector. Solid or crystallized Xenon detector can be ideal for dark 
matter searches, whereas a crystallized Xenon detector would be necessary to have channeling. At present, 
solid Xe crystals are used by the Solid Xenon R&D Project [lj . 

In this paper we compute the geometric channeling fractions of recoiling ions in solid Xe, Ar, and Ne 
crystals. Here "geometric channeling fraction" refers to assuming that the distribution of recoil directions 
is isotropic. In reality, in a dark matter direct detection experiment, the distribution of recoil directions 
is expected to be peaked in the direction of the average WIMP flow. At room temperature and pressure 
Xe, Ar and Ne are noble gases. At temperatures below 161.45 K, 83.80 K and 24.56 K respectively (at 
room pressure) they become solids. All of them form monatomic face-centered cubic (f.c.c.) crystals (see 
Appendix A). 

Our calculation is based on the classical analytic models developed in the 1960's and 70's, in particular 
by Lindhard and collaborators 21j . In these models the discrete series of binary collisions of the 

propagating ion with atoms is replaced by interactions between the ion and uniformly charged strings or 
planes. The screened atomic Thomas-Fermi potential is averaged over a direction parallel to a row or a 
plane of lattice atoms. This averaged potential is considered to be uniformly smeared along the row or plane 
of atoms, respectively, which is a good approximation if the propagating ion interacts with many lattice 
atoms in the row or plane by a correlated series of many consecutive glancing collisions with lattice atoms. 
We consider just one row or one plane, which simplifies the calculations and is correct except at the lowest 
energies we consider. 

We proceed in a similar manner as we did for Si and Ge in Ref. [10], and we refer to that paper for 
formulas and other details. There are several good analytic approximations of the screened potential. As in 
Ref [loj ]. here we use Moliere's approximation for the atomic potential, following the work of Hobler 21 1 and 



Morgan and Van Vliet 
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17j . In Moliere's approximation [221 ] the axial continuum potential, as a function 



of the transverse distance r to the string, is 

U(r) = (2Z 1 Z 2 e 2 /d) f(r/a) = E^f(r/a), (1) 

where E is the energy of the propagating particle and — 2Z\Zie 2 / '(Ed). Here Z\, Z-i are the atomic 
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numbers of the recoiling and lattice nuclei respectively, d is the spacing between atoms in the row, a is 
the Thomas-Fermi screening distance, a = 0.4685A(Z 1 1/2 + Z 2 1/2 )" 2/3 22, 3 and E = Mv 2 /2 is the 



22| for the continuum potential is 



kinetic energy of the propagating ion. Moliere's screening function 
/(£) = Si=i a iKo(Pi€)- Here Kq is the zero-order modified Bessel function of the second kind, and the 
dimensionless coefficients a t and fa are cti = {0.1,0.55,0.35} and fa — {6.0,1.2,0.3} [24j, for i = 1,2,3. 
The string of crystal atoms is at r = 0. In our case, E is the recoil energy imparted to the ion in a collision 
with a WIMP, 

and q is the recoil momentum. 

The continuum planar potential in Moliere's approximation [2^ . as a function of the distance a; perpen- 
dicular to the plane, is 

U p (x) = {2irnZ 1 Z 2 e 2 a) f p (x/a) = Eiplf p (x/a), (3) 

where ip 2 = 2irnZiZ2e 2 a/ E. Here n = Nd pc h is the average number of atoms per unit area, where TV 
is the atomic density and d pc h is the width of the planar channel, i.e. the interplanar spacing (thus, the 



average distance of atoms within a plane is d p = 1/ \/Nd pc h). The subscript p denotes "planar" and 
fp(0 — Y^=i ( a i/fa) ex P( — faO> where the coefficients cti and fa are the same as above. The plane is at 
x = 0. Also, the axial channel width which we call o? ac h is defined in terms of the interatomic distance d in 
the corresponding row as <i ac h = 1 /V Nd. Examples of axial and planar continuum potentials, generically 
called U, for Xe, Ar, and Ne ions propagating in the <100> axial and {100} planar channels of a Xc, Ar, 
and Ne crystal respectively are shown in Fig. [1] 

Lindhard proved that for channeled particles the longitudinal component of the velocity, i.e. the com- 
ponent along the direction of the row or plane of the velocity, may be treated as constant (if energy loss 
processes are neglected). Then, in the continuum model, the trajectory of the ions can be completely de- 
scribed in terms of the transverse direction, perpendicular to the row or plane considered. For small angle 
4> between the ion's trajectory and the atomic row (or plane) in the direction perpendicular to the row (or 
plane), the so called "transverse energy" 

E± = £sin 2 4> + U (4) 

is conserved. In Eq. 2] relativistic corrections are neglected. 

The conservation of the transverse energy provides a definition of the minimum distance of approach to 
the string, r m ; n (or to the plane of atoms x m - m ), at which the trajectory of the ion makes a zero angle with 
the string (or plane), and also of the angle ip at which the ion exits from the string (or plane), i.e. far away 
from it where U ~ 0. In reality the furthest position from a string or plane of atoms is the middle of the 
channel (whose width we call d ac h for an axial channel and d pc h for a planar channel). Thus, for an axial 
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Figure 1: Continuum axial (solid lines) and planar (dashed lines) potentials as function of the distance from the row or plane 
of lattice atoms, respectively, for (a) Xe ions, (b) Ar ions, and (c) Ne ions, propagating in the <100> axial and {100} planar 
channels of a Xe, Ar, and Ne crystal respectively. The screening radius shown as a vertical line is ax c Xc = 0.007808 nm for 
Xe, aiArAr = 0.01126 nm for Ar and a^eNc = 0.01370 nm for Ne (see App. A). 
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channel 

E ± = U{r min ) = EiP 2 + U{d ach /2). (5) 

The continuum model implies that the net deflection due to the succession of impulses from the peaks 
of the potential is identical to the deflection due to a force —U'. This is only so if the ion never approaches 
so closely any individual atom that it suffers a large-angle collision. Lindhard proved that for a string of 
atoms this is so only if U"(r) < 8E/d 2 , where the double prime denotes the second derivative with respect 
to r. Replacing the inequality in this equation by an equality defines an energy dependent critical distance 
r c such that r > r c for the continuum model to be valid. Morgan and Van Vliet [16| use 5 instead of 8 in this 
equation, because thi s ag rees better with their simulations of channeling in copper crystals. As in Ref. [lol ] 
and following Hobler [2l| (who studied channeling in Si crystals to be able to avoid it in ion-implantation) 
we use Morgan and Van Vliet's equation to define r c , i.e. 

U"(r c ) = ^E. (6) 

For a "static lattice", that here means a perfect lattice in which all vibrations are neglected, the critical 

n 

distance of approach r c is given in Eq. 2.15 of Ref. [101 ] . We use an approximate analytic expression for r c 
obtained by fitting a degree nine polynomial in the parameter y/a = (ZiZ2e 2 d/a 2 E) 1 / 2 to the exact solution 
of Eq. [6] The expression for r c obtained in this way (Eq. 2.15 of Ref. [lOj) is valid for recoil energies E > 3 
keV for Xe, and for E > 1 keV for Ar and Ne. 

Since r c is the smallest possible minimum distance of approach to the string of a channeled propagating 
ion for a given energy E, i.e. r m ; n > r c , and the potential U(r) decreases monotonically with increasing 
r, then C/(r m i n ) < U(r c ). Using Eq. [5l this condition translates into an upper bound on ip, if) < il) c {E) 
where tpc(E) is the critical channeling angle for the particular axial channel, i.e. the maximum angle 
the propagating ion can make with the string far away from it (in the middle of the channel) if the ion is 
channeled. We proceed similarly for planar channels to get the critical angles ipv(E). However, the treatment 
of planar channels presents some complications. 

The breakdown of the continuum theory for a planar channel is more involved than for an axial channel 
because the atoms in the plane contributing to the scattering of the propagating ion are usually displaced 
laterally within the plane with respect to the ion's trajectory. Thus the moving ion does not encounter atoms 
at a fixed separation or at fixed impact parameter as is the case for a row. To find the static critical distance 
x c of a planar channel we follow the procedure of defining a "fictitious string" , introduced by Morgan and 
Van Vliet QQ- 

They reduced the problem of scattering from a plane of atoms to the scattering from an 
equivalent fictitious row contained in a strip of width 2R centered on the projection of the ion path onto the 
plane of atoms. Along the fictitious row, the characteristic distance between atoms needs to be estimated 
using data or simulations which are not available for Xe, Ar and Ne crystals. As explained in Ref. jloj |. 
this characteristic distance depends on the width 2R of the strip considered. For R, Morgan and Van Vliet 
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used the impact parameter in an ion-atom collision corresponding to a deflection angle of the order of "the 



break-through" angle yJU v {Q) / ' E. This is the minimum angle at which an ion of energy E must approach 
the plane from far away (so that the initial potential can be neglected) to overcome the potential barrier at 
the center of the plane at x = (namely, so that E±_ = U P (Q)). For small scattering angles.the deflection 
angle 6 is related to the impact parameter, in this case R, as (see e.g. Eq. 2.1' of Lindhard (|) 

2E5 = -dU'(R), (7) 

where U' is the derivative of the axial continuum potential, and Morgan and Van Vliet define R by taking 



S = y/U p (0)/E. 

As described in Ref. [lOj, we decided to keep the Morgan and Van Vliet definition for R, because 
Hobler [21| found that it is in quite good agreement with the binary collision simulations and data of B and P 
ions propagating in Si for energies of about 1 keV and above. We used an approximate analytical solution for 
R (Eq. 2.22 of Ref. [10]) obtained by fitting a degree five polynomial in (lny) where y = Z\ Zi e 2 /a ^/ EU V (0) , 
to the exact numerical solution of Eq. [7] The planar critical distance x c we obtained (Eq. 2.23 of Ref. [lOj ) 
is valid for E < 7 GeV for Xe, E < 160 MeV for Ar, and E < 20 MeV for Ne. These conditions and those 
after Eq. [5] provide the energy ranges for which our channeling fraction estimates are valid. Within its range 
of validity, the percentage error of the analytic approximation we used for x c is less than 9%. 

The critical distances of approach in a non-static lattice depend on the temperature, through the vibration 
of the atoms in the lattice (thermal expansion effects are negligible, as shown in Appendix B of Ref. 
We use the Debye model to account for the vibrations of the atoms in a crystal. The one dimensional rms 
vibration amplitude u\ of the atoms in a crystal in this model (see Eqs. 2.27 and 2.28 in Ref. [ljj) is plotted 
in Fig. [5] for Xe, Ar, and Ne crystals as function of the temperature T up to their respective melting points. 
The crystals in the Solid Xenon R&D Project experiment will be operating at 77.2 K and higher [251 ]. 

In principle there are modifications to the continuum potentials due to thermal effects, but we take into 
account thermal effects in the crystal through a modification of the critical distances found originally by 
Morgan and Van Vliet 16] and later by Hobler 21 j to provide good agreement with simulations and data. 
For axial channels it consists of taking the temperature corrected critical distance r c (T) to be r c (T) = 



y^r^(E) + [ciMi(T)] 2 . For planar channels the situation is more complicated, because some references give 
a linear and others a quadratic relation between x c (T) and u\. Following Hobler [21| we use an equation 



similar to that for axial channels, namely x c {T) = y/ 'x^(E) + [c2U\(T)] 2 . The dimensionless factors c\ and C2 
were found to be numbers between 1 and 2 in the prior literature, for different crystals and propagating ions 
(see Ref. Q] and loj). Using the temperature corrected critical distances of approach r c (T) and x c (T) or 
the static lattice critical distances r c and x c (Eqs. 2.15 and 2.23 of Ref. [l(3])> we obtain the corresponding 
temperature corrected critical axial and planar channeling angles, ip c and ip? respectively. As c\ and C2 
increase, r c (T) and x c (T) increase and the corresponding critical angles decrease, respectively, at a given 
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Figure 2: Plot of u\(T) for Ne (thick solid line), Ar (thick dashed line), and Xe (thick dotted line) as a function of the crystal 
temperature. The plots are cut at the respective melting temperatures. 

temperature. Taking c\ = c% = corresponds to considering only a static lattice, with no temperature effects 
included. Examples of critical distances and angles are shown in Figs. |3]to[5j for c\ — C2 — c and c = 1 or 
c = 2 as indicated. Fig. [3] clearly shows how the critical distances and angles change with temperature for 
Xe ions propagating in the <100> axial channel of a Xe crystal, with temperature effects computed with 
c\ = C2 = 1. At small energies the static critical distance of approach is much larger than the vibration 
amplitude u\, so temperature corrections are not important. As the energy increases, the static critical 
distance of approach decreases, and when it becomes negligible with respect to the vibration amplitude m, 
the temperature corrected critical distance r c becomes equal to c\U\. Fig.0]shows the same as Fig.[3]but for 
c\ = C2 = 2. Figs. Efa) and|5^b) show the critical distances at several temperatures for Xe ions propagating 
in the {100} planar channel with c\ = C2 = 1 and c\ = C2 = 2 respectively. 

2. Channeling fractions 

In our model, a recoiling ion is channeled if the collision ion- WIMP happens at a distance large enough 
from the string or plane to which the ion belongs. Namely, channeling happens if the initial position of the 
recoiling motion is r, > r^min or Xi > a^mm f° r an axial or planar channel respectively. We define E± in 
terms of the initial recoil energy E of the propagating ion, the angle of the initial recoil momentum with 
respect to the particular string or plane of atoms <pi , and the initial position r,; . For an axial channel, the 
condition 

E ± (E,<t> i ,r i ) = U(r min ) <U(r c ), (8) 
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Figure 3: Static (thin dot-dashed line) and temperature corrected with c\ = C2 = c = 1 (thick lines) (a) critical distances of 
approach (and ui(T)) and (b) the corresponding critical channeling angles, for T = 77.2 K, T = 130 K, and T = 160 K as a 
function of the energy of propagating Xe ions in the <100> axial channel of a Xe crystal. 
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Figure 5: Static (thin dot-dashed line) and temperature corrected (thick lines) with (a) c\ = C2 = c = 1 and (b) c\ = C2 = c = 2 
critical distances of approach for T = 77.2 K, T = 130 K, and T = 160 K as a function of the energy of propagating Xe ions 
in the {100} planar channel of a Xe crystal. 



i.e. E±(E, <pi, r^ m i n ) = U(r c (E)) defines the distance r^ m i n in terms of r c . The equivalent formulas for axial 
channels define Xi jm - m in terms of x c . In Ref. ,5||, we derived an analytic expression for r^ m i n and a^min using 
Lindhard's approximation to the potential (see Eqs. 5.9 and 5.11 of Ref. [5j). 

We take the initial distance distribution (in rori respectively) of the colliding atom to be a Gaussian 
with a one dimensional dispersion U\, and we obtain the probability of channeling for each individual channel 
by integrating the Gaussian between the minimum initial distance and infinity (a good approximation to the 
radius of the channel; see Ref. [5j for details). The dependence of these probabilities on the critical distances 
enter in the argument of an exponential or an erfc function. Thus any uncertainty in our modeling of the 
critical distances becomes exponentially enhanced in the channeling fraction. This is the major difficulty of 
the analytical approach we are following. 

In order to obtain the total geometric channeling fraction we need to sum over all the individual channels 
we consider. The integral over initial directions is computed using HEALPix 26 j (see Appendix B of Ref. Q). 
Taking only the channels with lowest crystallographic indices, 100, 110 and 111, we have a total of 26 axial 
and planar channels (see Appendix A). We treat channeling along different channels as independent events 
when computing the probability that an ion enters any of the available channels. In Appendix D of Ref. [f| 
we showed that we can obtain an upper limit to the channeling probability of overlapping channels by 
replacing the intersection of the complements of the integration regions in Eqs. 5.2 and 5.4 of Ref. 5] with 
the inscribed cylinder of radius ?"min equal to the minimum of the ri jm - m or x^ m i n among the overlapping 
channels. We find that this method gives results practically indistinguishable from those obtained assuming 
that channeling along different channels are independent events. 

Fig. [S] shows the channeling probability computed for each initial recoil direction q plotted on a sphere 
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Figure 6: (color online) Geometric channeling probabilities as function of the initial recoil direction for (a) a 1 MeV Xe ion 
propagating in a Xe crystal at 160 K and (b) a 20 keV Ne ion propagating in a Ne crystal at 23 K (with c\ = C2 = 1}- The 
probability is computed for each direction and plotted on a sphere using the HEALPix pixelization. The colors/grayscale shown 
in the horizontal bar between black and white (blue and red online) corresponds to values between and 1 in increments of 
0.1. 

using the HEALPix pixelization for (a) a 1 MeV Xe ion propagating in a Xe crystal at 160 K and (b) a 20 
keV Ne ion propagating in a Ne crystal at 23 K (ci = c% = 1 is assumed for the temperature effects). The 
white and black (red and blue online) indicate a channeling probability of 1 and zero, respectively (see the 
colors/grayscale in the figure). 

Fig. [7] shows the channeling fractions for several individual channels of (a) Xe ions propagating in a Xe 
crystal at T — 77.2 K, (b) Ar ions propagating in an Ar crystal at T = 40 mK, and (c) Ne ions propagating 
in a Ne crystal at T — 40 mK (again with c\ — C2 = 1). The curves correspond to single axial or planar 
channels. Notice that at low energies channeling is dominated by axial channels, and at higher energies 
planar channels dominate. 

The geometric total channeling fractions of Xe, Ar, and Ne ions as function of the recoil energy are 
shown in Figs. [51 El and ITU1 respectively. For each crystal, we include three possibilities for thermal effects: 
(a) ci = C2 = 1, a reasonable middle ground, (b) maximum effects, i.e. c\ — C2 = 2, and (c) c\ — ci = 0, 
which corresponds to the unrealistic case of not having thermal effects in the lattice. 

Notice that we have not considered the possibility of dechanneling of initially channeled ions due to 
imperfections in the crystal. Any mechanism of dechanneling will decrease the fractions obtained here. 

As we see in Fig. [5] to [TU] the channeling fraction increases with energy, reaches a maximum at a certain 
energy, then has a dip and finally raises again. The maximum reflects the shape of the single channeling 
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Figure 7: Channeling fraction of (a) Xe recoils at T = 77.2 K, (b) Ar recoils at T = 40 mK, and (c) Ne recoils at T = 40 mK 
for single planar and axial channels, as a function of the recoil energy E in the approximation of c\ = C2 = 1. 



fractions, which all have maxima. These happen because the critical distances decrease with the ion energy 
E, which makes channeling more probable, while the critical angles also decrease with E, which makes 
channeling less probable. At low E the critical distance effect dominates, and at large E the critical angle 
effect dominates. At the maximum of the channeling fraction the axial channels dominate, in particular the 
channels [110] and [100] (as seen in Fig [7]). The dip and the raise result from having multiple axial and 
planar channels contributing to the channeling fraction. At lower E axial channels dominate and at higher 
E planar channels dominate. The dip happens at the cross-over between both types of channels, when as 
the energy increases the contribution of axial channels dies out and that of planar channels is increasing. 
This increase causes the subsequent raise in the channeling fraction as the energy increases further. 

As shown in Figs. [5] (a), |H] (a), andfTUl (a), the channeling fractions are never larger than 1% for Xe and 
Ar and never larger than 2% for Ne (with ci = c% = 1). 

Temperature effects computed with c\ = C2 = 1 are in the low range of what is found in other materials 
and propagating ions. Considering larger values of c\ and C2, the channeling fractions are smaller. Figs. [8] 
(b), IH1 (b), and [TU1 (b) show that with c\ — c-i — 2, the maximum channeling fraction for Xe ions at 160 
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Figure 8: Channeling fraction of Xe recoils as a function of the recoil energy E for T = 160 K (solid line), 130 K (dashed line), 
and 77.2 K (dotted line) in the approximation of (a) ci = C2 = 1, (b) c\ = C2 = 2 and (c) static lattice with c\ = c^ = 0. 



K, Ar ions at 83 K, and Ne ions at 23 K would be below 0.2%. However, since we do not know which are 
the correct values of the crucial parameters c\ and Ci for Xe, Ar, and Ne, we could ask ourselves how the 
channeling fractions would change if the values of these parameters would be smaller than 1 (although for 
other materials and propagating ions the values of these parameters were found to be between 1 and 2). 
The values of c\ and C2 cannot be smaller than zero, thus Figs. |8] (c), [9] (c), and ITOl (c) show the upper 
bounds on the geometric channeling fraction, obtained by setting c\ = C2 = 0, namely by neglecting thermal 
vibrations of the lattice (which make the channeling fractions smaller as T increases) but including the 
thermal vibrations of the nucleus that is going to recoil (which make the channeling fraction larger as T 
increases). Although it is physically inconsistent to take only the temperature effects on the initial position 
but not on the lattice, we do it here because using c± = C2 = 0, namely a static lattice, provides an upper 
bound on the channeling probability with respect to that obtained using any other non-zero value of c\ or 
C2- Even in this case, the channeling fractions cannot be larger than 5%. 

The critical angles depend on the temperature through the minimum distances of approach: as these 
increase with increasing temperatures, the critical angles decrease, what makes the channeling fraction 
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Figure 9: Channeling fraction of Ar recoils as a function of the recoil energy E for T = 83 K (solid line), 77.2 K (dashed line), 
and 40 mK (dotted line) in the approximation of (a) ci = C2 = 1, (b) c\ = C2 = 2 and (c) static lattice with c\ = C2 = 0. 



smaller. However, there is a second temperature effect which makes the channeling fractions larger as the 
temperature increases: the vibrations of the atom which collides with the dark matter particle. Depending 
on which of the two competing effects is dominant, the channeling fraction may either increase or decrease as 
the temperature increases. Increasing the temperature of a crystal usually increases the fraction of channeled 
recoiling ions, but when the values of ci and C2 are large (i.e. close to 2) so the critical distances increase 
rapidly with the temperature, the opposite may happen (see Fig. |H](b)). 

To conclude, let us remark that the analytical approach used here can successfully describe qualitative 
features of the channeling and blocking effects, but should be complemented by data fitting of parameters 
and by simulations to obtain a good quantitative description too. Thus our results should in the last instance 
be checked by using some of the many sophisticated Monte Carlo simulation programs implementing the 
binary collision approach or mixed approaches. 
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Ne ions, Ci=C2=l 



Ne ions, c\=C2=2 




Figure 10: Channeling fraction of Ne recoils as a function of the recoil energy E for T = 23 K (solid line) and 40 mK (dashed 
line) in the approximation of (a) c\ = C2 = 1, (b) c\ = C2 = 2 and (c) static lattice with c\ = C2 = 0. 
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Appendix A. Crystal structure and other data for Xe, Ar, and Ne 

Solid Xe, Ar and Ne have a face-centered cubic (f.c.c.) lattice structure with 4 atoms per unit cell. The 
lattice constant of Xe, Ar and Ne crystals are = 0.620 nm at T = 75 K 



27J, a(£ = 0.525 nm and 



a iat = 0-442 nm at T = 4.2 K and atmospheric pressure [28(. 

The atomic mass and atomic numbers of Xe, Ar, and Ne are Afxc = 131.29 amu, Mai = 39.948 amu, 
M Nc = 20.1797 amu, Z Xc = 54, Z Al = 18 and Z Nc = 10. 

The Thomas-Fermi screening distance for an ion expelled from a lattice site in the crystal scattering 
on another atom in the same crystal is axoXc = 0.4685(Z X { 2 + ZxeV 2/3 = 0.007808 nm for Xe, OArAr — 
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0.4685(Za{ 2 + Zi{ 2 )" 2/3 = 0.01126 nm for Ar and a NoN o = 0A685(Z^ C 2 + Z^ 2 )- 2 ? 3 = 0.01370 nm for Ne. 



To compute the interatomic spacing d in axial directions and the interplanar spacing d pc h ( "pch" stands 
for "planar channel") in planar directions, we have to multiply the lattice constant by the following cocfE- 



• Axis: <100>: 1 , <110>: 1/V2 , <111>: V3 

• Plane: {100} : 1/2 , {110} : 1/2^2 , {111} : l/VS 

The Debye temperatures of Xe, Ar and Ne are 9 Xc = 55 K, 9 Ar = 85 K and 9 Nc = 63 K Q, 
and the crystals in the Solid Xe R&D Project experiment will be operating at a temperature of 77.2 K or 



References 

[1] E. M. Drobyshevski, Mod. Phys. Lett. A 23, 3077 (2008) |arXiv:0706.3095l [phvsicB.ins-det]]. 
[2] R. Bernabei et ai, Eur. Phys. J. C 53, 205 (2008) |arXiv:0710!0288l [astro-phll. 
[3] R. Bernabei et ai, Eur. Phys. J. C 56, 333 (2008) | arXiv:0804^74TV l [astro-ph]]. 

[4] C. Savage, G. Gelmini, P. Gondolo and K. Freese, JCAP 0904, 010 (2009) |arXiv:0808!3607l [astro-phll. 
[5] N. Bozorgnia, G. B. Gelmini and P. Gondolo, JCAP 11, 019 (2010) |arXiv:1006!3Ti0l [astro-ph. CO] ]. 
[6] J. Lindhard, Kongel. Dan. Vidensk. Selsk., Mat.-Fys. Medd. 34 No. 14 (1965). 

[7] C. Savage, G. Gelmini, P. Gondolo and K. Freese, Phys. Rev. D 83 055002 (2011) |arXiv:1006!0972V 2 [astro-ph.CO]] . 

[8] F. T. Avignone, R. J. Creswick and S. Nussinov, arXiv:0807.3758 [hep-ph]; R. J. Creswick, S. Nussinov and F. T. Avignone, 
larXiv:1007!Q2T4V 2 [astro-ph.IM]. 

[9] N. Bozorgnia, G. B. Gelmini and P. Gondolo. larXiv:1101.28"76l [astro-ph.COl. 
[10] N. Bozorgnia, G. B. Gelmini and P. Gondolo, JCAP 11, 028 (2010) arXiv:1008.3676l [astro-ph.COll. 
[11] N. Bozorgnia, G. B. Gelmini and P. Gondolo, JCAP 11, 029 (2010) arXiv:1009.3325l [astro-ph.COll. 
[12] D. N. Balakishiyeva, R. Mahapatra, T. Saab and J. Yoo AIP Conf. Proc. 1274, 138 (2010). 
[13] J. Yoo, http://www.physics.ucla.edu/hep/dmlO/talks/yooxenon.pdf 

[14] Chapter 2 of "Ion Implantation", by Geoffrey Dearnaley, Amsterdam, North-Holland Pub. Co.; New York, American 
Elsevier (1973). 

[15] J. U. Andersen, Kongel. Dan. Vidensk. Selsk., Mat.-Fys. Medd. 36, No. 7 (1967). 

[16] D. V. Morgan and D. Van Vliet, Can. J. Phys. 46, 503 (1968); D. V. Morgan and D. Van Vliet, Radiat. Effects and 

Defects in Solids, 8, 51 (1971). 
[17] D. Van Vliet, "Channeling", D.V. Morgan ed., Wiley, London, p. 37 (1973). 
[18] J.U. Andersen and L.C. Feldman, Phys. Rev. B 1, 2063 (1970). 
[19] K. Komaki and F. Fujimoto, Phys. Stat. Sol. (a) 2, 875 (1970). 

[20] B. R. Appleton and G. Foti, "Channeling" in Ion Beam Handbook for Material Analysis, edited by J. W. Mayer and E. 

Rimini (Academic, New York), p. 67 (1977). 
[21] G. Hobler, Radiation effects and defects in solids 139, 21 (1996); G. Hobler, Nucl. Instrum. Methods Phys. Research 

(NIM) B 115, 323 (1996). 
[22] D. S. Gemmell, Rev. Mod. Phys. 46, 129 (1974). 
[23] J. H. Barrett, Phys. Rev. B 3 1527 (1971). 



cients 22j: 



higher 




16 



[24] G. N. Watson, "Theory of Bessel Functions" Cambridge U. P., Cambridge, England (1958). 

[25] J. Yoo, private communication. 

[26] K. M. Gorski et al, ApJ 622, 759 (2005). 

[27] D.R. Scars and H.P. Klug, J. Chem. Phys. 37, 3002 (1962). 

[28] D. G. Henshaw, Phys. Rev. Ill, 1470 (1958). 

[29] H. Fenichel and B. Serin, Phys. Rev. 142, 490 (1966). 

[30] E. R. Dobbs and G. O. Jones, Rep. Prog. Phys. 20, 516 (1957). 



17 



